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Abstract 

We derive a rigorous dynamical relation on aging phenomena — the aging 
relation — for Ising spin glasses using the method of gauge transformation. 
The waiting-time dependence of the auto-correlation function in the zero- 
field-cooling process is equivalent with that in the field-quenching process. 
There is no aging on the Nishimori line; this reveals arguments for dynamical 
properties of the Griffiths phase and the mixed phase. The present method 
can be applied to other gauge-symmetric models such as the XY gauge glass. 
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Slow dynamics is a key concept to study complex systems such as spin glasses, struc- 
ture glasses, polymers, superconductors, neural networks, and so on. It is one of peculiar 
properties to characterize the spin glass (SG) phase pHl. In the mean field theory, the non- 
ergodicity below the freezing temperature is significant to understand the nature of spin 
glasses @||. The aging phenomenon |J is a typical realization of slow dynamics especially 
in the SG |7HT0f. It was first observed for a metallic SG material CuMn in the zero-field- 
cooling (ZFC) process J7J. The relaxation of the isothermal remanent magnetization depends 
on the waiting time in which the sample is kept at constant temperature prior to the field 
application. Similar waiting time dependence was observed for other metallic SG materials 
and short-range SG materials (e.g. AgMn, AuFe, Fe^Mn^ .^TiOs, CdCr 3! Iii,f_ a ;S4, ■ ■ •), and 
on the thermo-remanent magnetization in the field-cooling (FC) process |8|-|10fl. 

Several attempts have been made theoretically to explain the aging by phenomenological 
arguments |lT|--[l3|. Cugliandolo and Kurchan [|TJ[] investigated the aging for the Sherrington- 



Kirkpatrick (SK) model analytically. They examined the auto-correlation function for the 



nonequilibrium process from a random state (ZFC from T = oo). Rieger investigated the 
aging in the 3D ±J Ising model by the Monte Carlo simulation. He measured the waiting- 
time dependence of the auto-correlation function from the all-up state (field quenching; see 
below) . 

In the present article, we derive a rigorous dynamical relation on aging phenomena for 
Ising SG models between two nonequilibrium processes shown in Fig. 1. The auto-correlation 
functions with a waiting time t w for these processes have a relation, 

+fw)>*] c = [(W(t+aJ'] c , (!) 

where K = J/k^T is the inverse temperature of the heat bath and K p gives the effective 
(inverse) temperature characterizing the randomness (see later). The process I is relating to 
the ZFC. At initial time t = 0, the system is kept in the equilibrium state of a temperature 
K p with zero field; hereafter, we use a terminology "temperature" for K and K p instead of 
"inverse temperature". The temperature is immediately changed (usually quenched) and 
the system relaxes in a heat bath with another temperature K in t > 0. The average for 
dynamical ensembles in this process is denoted by (• • -)^ p . The process II is relating to the 
field quenching (FQ) P,p~5f. The system starts from the all-up state F = (+, +,■••,+) at 



t — and relaxes in the same heat bath as in the process I in t > 0; the average is denoted by 
(• • •)#•. Since the all-up state provides the state in the strong-field limit, this represents the 
process with the field quenched from oo to zero at t = 0. Note that the FQ is not equivalent 
to the FC in which the applied field is more weak and is quenched after the waiting time. 
In the following, we show the derivation of Eq. ([!]) which we call the aging relation using 
the method of gauge transformation, |T^[19| and discuss the physical meaning of it. 



Since randomness and frustration make it difficult to examine SG systems analytically 
as well as numerically, only few things have been confirmed definitely. The method of gauge 
transformation is a powerful technique to derive exact results in the ±J or the Gaussian 
Ising spin glasses irrespective of the dimensionality. It provides the internal energy and an 
upper bound on the specific heat as non-singular functions of the temperature on a special 
line in the randomness-temperature phase diagram. This line is called the Nishimori line 
ITS| ] . Further, it provides a plausible argument for the absence of re-entrant transition from 
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the FM phase to non-FM one (SG in 3d) ||17|| . The method has been generalized to other 
random systems with various symmetries such as the XY gauge glass 
extended to treat dynamical systems . 
The Hamiltonian we consider is 



18 1 . Recently, it is 



U = JH(S;uj) 



-J LOijSiSj, 

(a) 



(2) 



where Si takes ±1, and S = (Si, 5*2, ■ ■ • , Sn) represents a configuration of the total N spins. 
The set uj = (uju, ■ • ■) represents a configuration of the total iVe bonds, and the summation 
is taken over all bonds; while we make no restrictions on the type or the dimension of the 
lattice, one may suppose usual nearest-neighbor interactions on the d- dimensional hypercubic 
lattice. For a particular bond configuration u>, the equilibrium distribution is defined by 
p eq (S; K, uj) = exp { — KTt(S; uj)}/Z(K, u>). The exchange interaction Jij = Jujij is a 
random variable. We treat both the ± J distribution and the Gaussian distribution with the 
variance unity. The average for bond configurations is denoted by [• • -] c . The general form 



of the bond distribution in a gauge-symmetric model is expressed as 



P(u>\K n 



DM 

Y(K P 



exp{-i^(F;u,)}, 



(3) 



where F = (+,+,•■■, +). Functions u^, K p , D(u>) and Y(K p ) are summarized in table 
|. In both distributions, K p controls the randomness; K p = and oo correspond to the 
most random case and the non-random case, respectively. The Nishimori line is located on 
K = K p . 

Since the Ising system has no intrinsic dynamics, we consider a Markov process for each 
bond configuration: The density of state obeys the master equation [EDI, 



-jiMS) = £r(s|sV(s'). 

The solution of the master equation is formally given by 

p,(S) = ]T <S|e* r |S'> po(S'). 



(4) 



(5) 



The matrix element (S|e tr |S') plays a role of a Green's function for a time interval t. The 
matrix T is composed of non-negative off-diagonal elements, and satisfies the detailed balance 



and the conservation of the probability. We consider both the Metropolis [21| and the 



Glauber dynamics ||22|| . The detailed expressions of T for these dynamics are shown in Ref. 
0. 



The gauge transformations for functions of S and uj are defined by 



S; 



SiCTi, 
&ijO'i<Jj, 



(6) 
(7) 



where a = (a\ , a<i , ■ ■ ■ , a^) is an arbitrary state of iV Ising spins. The Hamiltonian is 



invariant under the transformation Ua-Va-, 
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U^H^u) = H{S;u). (8) 

Another important property in the static case is the invariance of summation (or integral) 
for S and u>; 

£••• = E^---> (9) 

s s 

Note that we use the terminology "gauge invariant" only for functions of the set uj invariant 



under V^. The bond average of a gauge- invariant function Q(us) can be expressed as p8|J19 



To treat dynamical systems, another transformation for S' is necessary; 

V' a : fljj — ► flfo. (12) 

We have shown the invariance of time evolution 

VJJ' V V„ (S|e tr |S'> = (S|e tr |S'>, (13) 

for the Metropolis and the Glauber dynamics |19|| . This leads us to the relations obtained 
previously; 

rnrni = [(SiWswfci, (M) 

[W))*] c = [<H(t)) K K % (15) 

Equation (TJ) is a generalization of the fluctuation-dissipation theorem in the region far 
from equilibrium. 

We examine the aging by analyzing the waiting-time dependence of nonequilibrium auto- 
correlation functions |f[H|l5] for the processes I and II instead of the remanent magnetization 
like in experiments. Since the auto-correlation function is a fundamental quantity to exhibit 
the dynamical structure in the system, we believe it detects a typical aging feature of the 
process if any. In the language of the master equation, they are expressed as 

+ *w))* 

= E ^(S 2 |e* r |S 1 )5 li (S 1 |e^ r |F), (16) 
S1S2 

+ tw)>X 

= E^(S 2 |e' r |S 1 )5 H (S 1 |e'- r |S )p cq (S ;/^,a;). 
S0S1S2 

(17) 

From the invariances (|3l)-(|i~0|) and (|T3"D, it is easily seen that (Si(t w ) Si(t+t w ))^ is transformed 

as 
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V^S^S^t + t w ))l = (SiitJSiit + t w ))£, (18) 



while (Si(t w )Si(t + t w )) K p is gauge invariant; (• • expresses the average for the process 
starting from a fixed state cr. Note that the probability distribution (H) is transformed as 

V„P(u;K p ) = J( ">f^ M) p^;K p ,u,). (19) 



Using Eqs. (PD , (p~8| ) and (0), we derive the relation ([]]) for any waiting time t w , any time 
interval t, any temperature K and any degree of randomness K p . We call it the "aging 
relation" , since it relates aging phenomena in two distinct processes whatever waiting-time 
dependence is essential for the aging. Equation flU]) is a special case (t w = 0) of it. 

The aging relation contains parameters K and K p characterizing the relaxation in t > 0. 
One can examine the physical meaning of it for each phase by choosing (K, K p ) appropri- 
ately. A typical phase diagram in randomness-temperature plane for Ising spin glasses is 
shown in Fig. 2; possible Griffiths phase and mixed phase are indicated. Note that the 
initial temperature of the process I is always located on K = K p . In the SG phase where 
K p indicates a high temperature (small K p ), the process I expresses the ZFC observed in 
experiments. Therefore, the aging relation relates aging phenomena for the ZFC and the 
FQ processes. The same relation was pointed out in real experiments P,l0|, in which the 



waiting-time dependence of the remanent magnetization in the ZFC process seemed equiv- 
alent with that in the FQ process. Since the aging relation reveals just the equivalence of 
the gauge-invariant dynamical structure, the amplitude of the magnetization which is not 
gauge invariant does not have such relation, while dynamical behavior itself may have it 
as observed in these experiments. At a glance, the aging relation expresses a trivial fact, 
since both initial states are located in the PM phase and both systems relax into the SG 
phase. However, the equivalence at any waiting time, which means the equivalence of the 
dynamical structure at any stage of relaxation, is non-trivial. Further investigations in this 
direction would be helpful to understand the dynamical structure of the system. 

Next, let us consider the aging relation on K = K p (Nishimori line), where the temper- 
ature keeps a constant at any time in the process I. Then, the rhs of Eq. (Jl|) reveals the 
equilibrium auto-correlation function, which is independent of the waiting time t w : 

(S l (t w )S l (t + t w ))l] c = [<£(0)Si(f))£] c . (20) 

(• • denotes the dynamical average for the equilibrium process. Equation (|20|) is derived 
from the fact that the equilibrium distribution p eq (K) is an eigenstate of T with zero eigen- 
value. Therefore, the auto-correlation function in the FQ process is independent of the 
waiting time suggesting the absence of aging on the Nishimori line. Further, it is noted 
that Eq. (^) provides efficient Monte Carlo calculations for equilibrium auto-correlation 
functions; the equilibration can be omitted in the FQ process. 

The aging has been considered an inherent property in the SG phase from experimental 
f8HT0| and theoretical [|lTl- [l3|Jl5[| view points. Since the aging phenomenon is a typical 



realization of the complex phase space for slow dynamics, it could occur in other complex 
phases. In such a phase, the aging would also be inherent, which means that it occurs whole 
in the phase when it is observed in some part of the phase. It has been pointed out that there 
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exists a dynamically singular phase called the Griffiths phase P,|23|-p6| between the critical 
temperature of the pure system and the phase boundaries of low temperature phases (the 
FM or the SG) — see Fig. 2. Concerning to the Griffiths phase, two cases can be considered 
to understand the above result if the aging is an inherent property, (a) There is no aging 
whole in the Griffiths phase, if the Nishimori line intersects the Griffiths phase as in Fig. 2. 
Then, even if slow dynamics is observed in the Griffiths phase, it is quite different from that 
in the SG phase |27|] . However, the region of the Griffiths phase has not yet been determined 



definitely: Figure 2 was proposed just following an analogy with the dilute ferromagnet [28 



One can also consider that (b) there is no Griffiths phase at least around the Nishimori line. 
This allows both the existence of the Griffiths phase below K = K p and the absence of it. 
Another possibility is that the aging is not inherent at least in the Griffiths phase. Although 
it is not obvious which is correct in the present framework, the above result restricts the 
aging phenomenon and the region of the Griffiths phase in future investigations. 

Following the results of the SK model and experiments one should consider the 
possibility of the mixed phase between the FM and SG phases - see Fig. 2. It is reasonable 
to consider that the aging occurs in the mixed phase, since the SG feature in the mixed phase 
provides typical slow dynamics which reveals the aging. Let us consider the temperature 
below the multicritical point on K = K p , where the spontaneous magnetization appears. In 
such a region, it is not clear that the FQ process is appropriate for the observation of the 
aging feature, since the initial all-up state is not so different from the final equilibrium state 
with broken symmetry. As seen above, the amplitude of the magnetization is not important 
in the aging phenomenon. While the up-down symmetry is broken in the FQ process, 
the dynamical behavior is equivalent with that in the equilibrium (symmetric) process - 
Eq. ( flop . Thus, we assume that the FQ process provides an aging feature even in such 
a symmetry-broken region. If the aging is inherent in the mixed phase, the relation ( |2"DD 
indicates that the Nishimori line does not enter the mixed phase. This restricts the topology 
of the phase diagram, and is consistent with the result of the SK model. 

In summary, we derive the relation ([!]) on the aging phenomena for two nonequilibrium 
processes. This provides some restrictions to the Griffiths and the mixed phases. Similar 
relations can be derived for other gauge-symmetric quantities such as the SG susceptibility, 
and in other processes like in the simulated annealing. The present theory is applicable to 
various systems: Any gauge-symmetric distribution for randomness instead of the Gaussian 
and the ± J, any dynamics which satisfies Eq. (|T3|) instead of the Glauber and the Metropolis, 
any dimensionality and lattice (e.g. the SK model). It can be extended to other gauge- 
symmetric systems such as the XY gauge glass |fL8|l ; H = —JJ2(ij) cos(0j — <pj + Aij). It has 
been pointed out that granular systems of type-II superconductors in a magnetic field are 
expressed this model. There exists a similar relation |29| in such materials. 

The author would like to thank H. Rieger, E. Vincent, J. M. Hamman, I. A. Campbell, 
S. Miyashita, K. Nemoto, J. Kurchan and L. F. Cugliandolo for valuable discussions and 
comments. 
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FIG. 1. Illustrations of two nonequilibrium processes, I and II. Between t w and t + t w , the 
correlation is measured. 




TABLES 

TABLE I. The summary of variables and functions appeared in the bond distributions. Nb is 
the number of bonds, p is the concentration of + J bonds in the ± J distribution, and Jq is the 



center of distribution 


in the Gaussian distribution with the variance unity. 








±J 


Gaussian 


U>ij 




±1 


any real values 


K p 




iln(p/l -p) 


Jo 


Y(K P ) 




(2coshK p ) AfB 


exp (\N B K P 2 ) 


D(u) 




1 


exp ( - \ <4 3 ) 



9 



